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We introduce the concept of nonlinear graphene metasurfaces employing the controllable inter¬ 
action between a graphene layer and a planar metamaterial. Such hybrid metasurfaces support 
two types of subradiant resonant modes, asymmetric modes of structured metamaterial elements 
(” metamolecules”) and graphene plasmons exhibiting strong mutual coupling and avoided dispersion 
crossing. High tunability of graphene plasmons facilitates strong interaction between the subradiant 
modes, modifying the spectral position and lifetime of the associated Fano resonances. We demon¬ 
strate that strong resonant interaction, combined with the subwavelength localization of plasmons, 
leads to the enhanced nonlinear response and high efficiency of the second-harmonic generation. 
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Introduction. Unique electromagnetic properties of 
graphene are paving a way towards their applications 
in novel optical devices [1-9]. One of the main bene¬ 
fits of using graphene stems from the ability to control 
its optical response by changing the chemical potential 
via electrostatic or chemical doping [10-14]. While the 
electromagnetic response of bare graphene is considered 
to be strong, it is still not sufficient for a majority of 
practical applications. One of the promising approaches 
to enhance optical response of graphene is to place it 
into close proximity with metallic and dielectric photonic 
structures, and this approach has been employed to de¬ 
sign highly tunable infrared and THz metasurfaces en¬ 
abling both amplitude and phase modulation [15-21]. In 
addition, doped graphene exhibits a plasmonic response, 
and therefore it can be patterned to form tunable planar 
plasmonic metasurfaces operating at infrared and THz 
frequencies [9, 22, 23]. 

In this Letter, we discuss a novel type of Fano-resonant 
hybrid metasurfaces composed of a graphene layer and 
a planar metamaterial where two types of subradiant 
modes of different origin, asymmetric plasmonic modes 
of metallic metamolecules and graphene plasmons, are 
strongly coupled. We show that such metasurfaces 
support cascaded Fano resonances [24, 25] originating 
from the interference of subradiant plasmonic modes in 
graphene and localized modes of metamolecules. We 
demonstrate that such hybrid graphene metasurfaces 
provide strong tunability and dramatic field enhancement 
giving rise to enhanced nonlinear response. 

Model and parameters. We consider a metasur¬ 
face composed of a square lattice of metallic split-ring- 
resonators (SRRs) placed on top of a graphene layer, 
as shown in Figs. l(a,b). The electromagnetic response 
is calculated by using the finite-element method solver 
COMSOL Multiphysics with a graphene layer modelled 
as a surface current. At low frequencies hw < 2£p, 
graphene is described by a frequency-dependent surface 


conductivity [26] 
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where e = — \e\ is electron’s charge, £p = 
is the Fermi energy, n is the doping electron density, 
Vp ~ c/300 is the Fermi velocity, and Ti ntra = 0.15 ps 
is the relaxation time. At frequencies fuj < £p , the 
Drude-like conductivity originating in intra-band tran¬ 
sitions dominates and graphene supports highly confined 
extremely short-wavelength p-polarized plasmons, whose 
dispersion can be tuned by changing the Fermi level of 
graphene £p. Considering the normal incidence of a 
^/-polarized plane wave, we calculate numerically both 
transmission and reflection spectra for different values of 
the chemical potential. 

Without a graphene layer, the metasurface becomes 
a conventional two-dimensional metamaterial of double 
SRRs, supporting symmetric (superradiant) and asym¬ 
metric (subradiant) modes that correspond to the in- 
phase and out-of-phase current distribution along two 
arms of the metamolecule [see Fig. 2(a)]. When the 
symmetry of the metamolecules is reduced, the interfer¬ 
ence between these two modes results in a characteristic 
Fano-like asymmetric lineshape in far-field spectra, the 
transition illustrated by black dotted and dashed lines in 
Fig. 2(c) [27-30]. By introducing a graphene layer, we 
can excite graphene plasmons whose wavelength can be 
tuned to match one of the characteristic scales of the 
metasurface. Note, however, that even without plas¬ 
monic modes of graphene tuned to the spectral range 
of interest, the surface conductivity of graphene modifies 
electromagnetic interactions on the metasurface leading 
to the overall blue shift of the resonances [16]. 

Due to the subwavelength nature of graphene plas¬ 
mons, it is natural to match their wavelength to the 
smallest scale associated with the gaps of the split-ring 
resonators, which effectively become a cavity for plas¬ 
mons [18]. Here we are interested in graphene plasmons 
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Figure 1. (Color online) (a) Artistic view of the light in¬ 
teraction with a hybrid metasurface created by a lattice of 
asymmetric gold SRRs placed on top of a graphene layer, 
(b) Top view and size definitions of a unit cell of the hybrid 
metasurface. 


forming a subradiant mode, and such mode can be de¬ 
signed from two plasmons localized within the SRR gaps 
and interacting with each other through SRR arms. This 
interaction gives rise to a new asymmetric subradiant 
mode with the plasmons within two gaps being out of 
phase [see Fig. 2(b)]. Owing to their matching symme¬ 
tries, this additional subradiant mode can efficiently cou¬ 
ple to the asymmetric mode of SRRs when their frequen¬ 
cies are tuned [see Fig. 2(c)], resulting in a cascaded dou¬ 
ble Fano resonance in the far-field spectra [31-33]. More¬ 
over, since the wavelength of graphene plasmons can be 
altered by changing the Fermi level, this enables dynam¬ 
ically controllable interaction between the two subradi¬ 
ant modes making one of them less pronounced. Direct 
numerical calculations of the far-field spectra shown in 
Fig. 3(a) confirm that the hybrid metasurface supports 
one superradiant (broad peak) and two subradiant reso¬ 
nances associated with the Fano resonances, one of them 
being highly tunable. Profiles of two subradiant asym¬ 
metric modes are shown in Fig. 2(b) for the case of their 
strong hybridization. Changing the frequency of the in¬ 
cident plane wave allows to see clearly the features of 
two resonances in near-field (see Supplementary Mate¬ 
rial [34]). 

Analytical approach. The basic principle of the forma¬ 
tion of the cascaded Fano resonances can be illustrated 
by applying the coupled-mode theory [28]. In this ap¬ 
proach, each mode of the metasurface is described by a 
complex amplitude (S', A, and P), corresponding to the 
bright symmetric and dark asymmetric modes of SRR, 
and a dark asymmetric graphene plasmonic mode, re¬ 
spectively. All modes are normalized in such a way that 
\S\ 2 , |A| 2 , |P| 2 are the energy densities. The amplitudes 
of the incident P + and reflected E~ waves are also nor¬ 
malized such that |p( +,- )| 2 represent the incoming and 
outgoing power fluxes. Then, the system dynamics is 



Figure 2. (Color online) (a,b) Near-field images of the reso¬ 
nances shown by the E z field component without and with 
a graphene layer, respectively, (c) Numerically calculated re¬ 
flection spectrum at Ef — 0.307 eV (red circles) along with 
the fitting obtained by the coupled-mode theory (blue solid 
curve). For the cases of symmetric and asymmetric split- 
rings, the spectra without graphene are shown by black dotted 
and dashed lines, respectively. Structure parameters are: the 
substrate permittivity £2 — 1.35 2 , P x m 2.5 pm and P y = 4.5 
/im, the parameters of the gold SRR are L y — 3 pm, g — 310 
nm, W— 0.3 pm, W y — 1 pm, W\ — 1.28 pm, W 2 — 0.6 pm, 
and gold thickness is h = 40 nm. 


described by a set of coupled equations 

' S = (- iuj s - 7 s)S + igiA + ik s E ; + , 

< A = {—iuiA — 1 a)A + ig\S + ig^P , ( 2 ) 

. P = {~iup ~ 7 p)P + *52^4 , 

where and are bare frequencies and decay rates 
of the modes, gi^ characterize the interaction strengths, 
and Kg determines coupling of the symmetric mode to an 
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Figure 3. (Color online) (a) Numerically calculated transmis¬ 
sion spectra of the metasurface compared to (b) analytical 
results. Black dotted curves in panels (a,b) illustrate the an¬ 
ticrossing between the dark modes. 


incident plane wave. In general, the total decay rate is 
composed of two parts 7 s,a,p — Js^ap + Ys'a Ohmic 
(resistive) and radiative loss, respectively. It is assumed 
that only the symmetric mode is driven by the time- 
harmonic incident field, while the asymmetric modes can 
be excited only indirectly through coupling. Additional 
constraints should be imposed on the coefficients to sat¬ 
isfy the energy conservation and time-reversal: 7 ™p = 0 , 

M = E^ = E +, where the numeric subscript 

(1 or 2 ) refers to the media surrounding the metasur¬ 
face. For simplicity, the structure is assumed to be top- 
bottom symmetric, which neglects substrate effect justi¬ 
fied by small optical contrast (the reflection coefficient 
without resonances r = 0 ) so that E~^ = E^i) + 
Complex reflection can then be found from expression 
r = E^jEt = {ix s S)/E+ and assumes the form of a 
” nested” Fano resonance: 
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which clearly reflects cascaded coupling of the plasmonic 
asymmetric mode to the radiative continuum via the 
asymmetric and symmetric modes of SRRs. The trans¬ 
mission coefficient can then be found from the equation 
t = E^/Et = 1 + r. 

A change in the graphene chemical potential leads to 
variation of the spectral position of the plasmonic res¬ 
onance ujp. Due to strong near-field coupling of two 
subradiant modes their interaction reveal anticrossing of 
associated Fano resonances, shown in Fig. 3(a). The 
analytical model nicely reproduces this behavior. In 
Fig. 3(a) transmission calculated with the use of COM- 
SOL Multiphysics is superimposed with the dispersion 
curves found from the coupled-mode theory, with pa¬ 
rameters recovered from the fitting. The spectral po¬ 


sitions of the noninteracting modes are shown by dashed 
lines. Here the graphene plasmon branch is obtained 
from the standing-wave condition 27 t// c sp ~ g , where 
k sp satisfies the dispersion of the p-polarized plasmons 
on the graphene layer on top of dielectric substrate 

1 £ 2 ATTCrtuj) 

£ 2 : —~^=^= H- = = —% -. Interac- 

tion between the resonances leads to their hybridization 
with dispersion branches of the normal modes shown by 
dotted lines. 

Since the resonant frequency of graphene plasmons de¬ 
pends strongly on doping, the plasmonic mode intersects 
the subradiant SRR mode. Their interaction results in 
the characteristic avoided crossing behaviour with the 
frequency splitting between two modes. The results 
of the numerical and analytical calculations are plotted 
alongside in Fig. 3(a) and Fig. 3(b), respectively, where 
one can see that the tunable anti-crossing behavior de¬ 
scribed by these approaches agrees perfectly well. 

It is interesting to mention that, in a sharp contrast to 
conventional tunable plasmonic resonances in graphene, 
the subradiant nature of the graphene resonance consid¬ 
ered here allows not only tuning its spectral position, 
but also coupling to a radiative continuum. Indeed, as 
demonstrated in Fig. 3, the intensity of the asymmetric 
plasmonic mode fades away as one moves away from the 
SRR asymmetric mode, evidencing the reduction of the 
radiative coupling caused by the suppression of the cas¬ 
caded interaction between the modes. Therefore, the pro¬ 
posed metasurface offers an unprecedented control over 
its scattering characteristics, including spectral position, 
radiative coupling strength, and phase, unattainable in 
conventional designs of Fano-resonant metasurfaces. The 
ability to control radiative coupling can be of special in- 



Figure 4. (Color online) Numerically calculated absorption by 
the metasurface for graphene doping Sf = 0.27 eV (purple 
circles) and Sf = 0.307 eV (red circles). The insets show 
the field distributions in the SRR gap in the vicinity of two 
main absorption peaks associated with the excitation of Fano 
resonances. 
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terest for applications where the tunable near-field en¬ 
hancement is a key requirement. To demonstrate the 
possibility of a control over strength of light-matter inter¬ 
action, below we consider two examples: tunable absorp¬ 
tion enhancement and nonlinear second-harmonic gener¬ 
ation. 

Tunable absorption. Frequency-dependent absorption 
is plotted in Fig. 4 for two levels of graphene doping, with 
the insets showing the field distribution in the gap in the 
vicinity of two absorption peaks stemming from excita¬ 
tion of subradiant modes. The enhanced absorption due 
to excitation of hybrid graphene plasmon-quadrupole res¬ 
onances clearly correlates with near-field enhancement. 
As expected, for the case of poorly interacting subra¬ 
diant modes, corresponding to the case of lower Fermi 
energy (£p = 0.27 eV, shown by purple circles in Fig. 4), 
the asymmetric plasmonic mode has smaller amplitude 
and the loss are not as prominent as for the asymmetric 
SRR resonance. For higher doping level (£p = 0.307 eV, 
shown by red circles in Fig. 4) the two dark modes are 
strongly hybridized leading to a higher absorption rate 
by the tunable subradiant mode. This confirms that the 
system considered here can be promising for various ap¬ 
plications where tunable strong field enhancement is in 
demand. 

Nonlinear effects. Being a strongly nonlinear mate¬ 
rial, graphene is of a great interest for nanoscale non¬ 
linear optics [35], especially in combination with plas¬ 
monic effects [36-38]. Extreme surface confinement of 
electromagnetic fields due to excitation of plasmons is 
known to result in enhanced nonlinear response. Even 
stronger field confinement in hybrid metasurfaces pro¬ 
posed here can further boost the efficiency of various 
nonlinear optical effects [29, 39, 40]. In particular, the 
local field enhancement can be exploited for nonlinear 
frequency conversion. As the demonstration of enhanced 
nonlinear effects associated with the resonant excita¬ 
tion of subradiant plasmonic mode in graphene, here we 
study the second-harmonic generation (SHG) caused by 
a quadratic nonlocal nonlinearity of graphene. 

The nonlinear surface current induced in graphene by 
the incident wave of the fundamental frequency (FF) 
co>, at the double frequency 2cj has the form ji(2cj) = 
a ijki ( CJ ) Ef V k Ef , where i,j, k,l are the Cartesian com¬ 
ponents. It depends on the gradient of FF electric field 
at the graphene surface, and it can be expressed through 
the second-order conductivity tensor [41], 
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which is derived within the quasi-classical approach 
based on the Boltzmann equation [42-45], applicable at 
low frequencies Huj < £p corresponding to the plasmonic 
regime. In what follows, we neglect nonlinearity in metal. 

The nonlinear response of the metasurface is simulated 
in COMSOL Multiphysics with the use of two coupled 
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Figure 5. (Color online) (a) Nonlinear source and (b) SH 
fields generated by the hybrid metasurface; the z-component 
of the electric field is evaluated 40 nm above the graphene 
layer (left), and a side-view of the x-component evaluated in 
the middle of the SRR gap (right). 


electromagnetic models assuming undepleted pump field. 
As the first step, we perform simulations at the funda¬ 
mental (pump) wavelength 10.65 /im (matching one of 
the radiation lines of CO 2 laser) and stimulating excita¬ 
tion of graphene plasmon resonance in the SRR gap for 
the doping level £p = 0.31 eV. Next, the nonlinear cur¬ 
rent induced at every point on the graphene surface for 
the second-harmonic (SH) frequency is calculated. This 
current is employed as a nonlinear source for the electro¬ 
magnetic simulation and the SH fields generated by this 
current are obtained. 

As expected, SHG predominantly arises from the gap, 
where the field of the asymmetric plasmonic mode is lo¬ 
calized and strongly enhanced. The calculated nonlinear 
current generated by this subradiant mode is found to 
have the dominant x-component driving the entire meta¬ 
surface at the SH frequency. As confirmed by numer¬ 
ical simulations shown in Fig. 5, this results in a non¬ 
vanishing nonlinear dipole moment aligned along x axis 
within the SRR gaps, which leads to the generation of 
a x-polarized plane wave in the far-held region. It is in¬ 
teresting that the SH wave appears to be cross-polarized 
with respect to ^/-polarized incident plane wave at the 
fundamental frequency, indicating the nonlinear polar¬ 
ization conversion. We find that due to the strong held 
enhancement associated with two complimentary trap¬ 
ping mechanisms of electromagnetic held, namely the 
plasmonic conhnement and a subradiant character of the 
plasmonic mode, the structure exhibits high efficiency of 
nonlinear conversion ~ 10“ 6 at the incident wave inten¬ 
sity of 1 MW/cm 2 even without matching to any reso¬ 
nances at the SH frequency. The animation illustrating 
time dependence of the SH held component E x radiated 
by the metasurface at the frequency of graphene plas¬ 
monic Fano resonance in the plane perpendicular to the 
SRR gap can be found in Supplementary Material [34] . 

Conclusions. We have proposed a novel type of 
graphene-based metasurfaces exhibiting cascaded Fano 
resonances originating from the coupling of standing 
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graphene plasmons and the modes of structured meta¬ 
surfaces. We have shown that such hybrid metasurfaces 
allow high tunability of the subradiant modes facilitated 
by controllable intermodal coupling. Such tunable in¬ 
teraction between subradiant modes, achieved though 
graphene doping, enables controlling the spectral posi¬ 
tion and lifetime of the Fano resonances. By utilizing two 
complimentary mechanisms of light confinement origi¬ 
nating in the subwavelength localization and subradi¬ 
ant nature of the engineered graphene plasmonic modes, 
we have shown tunable and enhanced light-matter in¬ 
teraction. In particular, a significant enhancement of 
light absorption and nonlinear effects have been demon¬ 
strated. The proposed concept envisions tunable non¬ 
linear response of either graphene or any nonlinear ma¬ 
terials sputtered on top of the metasurfaces supporting 
cascaded Fano resonances. The demonstrated possibility 
to enhance and switch the second-harmonic generation 
by graphene gating can be used in various applications, 
e.g. for the efficient generation of short pulses of the 
second-harmonic field. 
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